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ABSTRACT 

We present a comparison between simulation results and X-ray observational data on the evo- 
lution of the metallicity of the intra-cluster medium (ICM). The simulations of galaxy clusters 
have been carried out using a version of the TreePM-SPH GADGET-2 code that includes a de- 
tailed model of chemical evolution, by assuming three different shapes for the stellar initial 
mass function (IMF). Besides the Salpeter (1955) IMF, we used also the IMF proposed by 
Kroupa (2001) and the top-heavier IMF by Arimoto & Yoshii (1987). We find that simula- 
tions predict significant radial gradients of the Iron abundance, Zp c , which extend over the 
whole cluster virialized region. Using the Salpeter IMF, the profiles of Zp have an amplitude 
which is in a reasonable agreement with Chandra observations within 0.2i?5oo- At larger radii, 
we do not detect any flattening of the metallicity profiles. 

As for the evolution of the ICM metal abundance out to z = 1, it turns out that the results 
based on the Salpeter IMF agree with observations. We find that the evolution of Zp e in 
simulations is determined by the combined action of ( i) the sinking of already enriched gas, 
( ii) the ongoing metal production in galaxies and ( Hi) the locking of ICM metals in newborn 
stars. As a result, rather than suppressing the metallicity evolution, stopping star formation at 
z = 1 has the effect of producing an even too fast evolution of the emission-weighted ICM 
metallicity, with too high values of Zp c at low redshift within 0.2i?2oo- Finally, we compare 
simulations with the observed rate of type-la supernovae per unit B-band luminosity (SnUs). 
We find that our simulated clusters do not reproduce the decreasing trend of SnUs at low 
redshift, unless star formation is truncated at z = 1. 

Key words: Cosmology: Theory - Methods: Numerical - AT -Rays: Galaxies: Clusters - 
Galaxies: Abundances - Galaxies: Intergalactic Medium 



1 INTRODUCTION 

The high quality X-ray observations of galaxy clusters from the 
current generation of X-ray satellites are allowing now to trace 
in detail the pattern of t he metal enrichment of the intra-clu ster 
medium of galaxies (e.g.,[Mushotzky 2004; Werner et al. 20(|| for 
reviews). In turn, this information is inextricably linked to the his- 
tory of formation and evolution of the g alaxy populatio n as ob- 
served in the optical/near-IR band (e.g., iRenzin i 2004 and ref- 
erences therein). A number of independent observations have es- 
tablished that significant radial gradients of the Iron abundance 
are present in the ce ntral regions, 0.1J?5qq, of relaxed clus- 
ters and groups (e.g..lbe Grandi et alj|2004 Ivikhlinin et~ai]|2005l ; 
iRasmussen & Pon man 2007|), with enhancement of the metallicity 
associated to the brightest cluster galaxies (BCGs), while no evi- 
dence has been found that these gradients extend at larger cluster- 



centric distances (e.g.. ISnowden et al]|2008l) . Furthermore, deep 
exposures with the Chandra and XMM-Newton satellites have 
now opened the possibility of tracing the evolution of the ICM 
metal content within the central regions out to the larg est redshifts, 
z ~ 1 .3, where clusters have bee n identified so far. lBalestra et al.l 
(2007) and Maughan et al. (2008) have analysed fairly large sam- 
ples of distant clusters, extracted from the Chandra archive and 
found that the metallicity of the ICM within the central cluster re- 
gions has increased by about 50 per cent since z ~ 1. 

This positive evolution of the ICM metallicity in the central 
cluster regions is apparently in contradictio n with the lack of sig - 
nificant star formation at low redshift ( e.g., iRaffertY et al]|200d) . 
Based on a phenomenological approach. lEttoril 120051) showed that 
the evolution of the ICM metallicity is in line with the expectations 
from the observed cosmic rates of supernova (Sn) explosions and of 
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star formation. iLoewensteinl d2006h combined observations of the 
evolution of the ICM metallicity with data on the Sn rates and star 
formation rates to infer the relative role played by type la and II Sn 
(Snla and SnII hereafter). 

A different approach was pursued by other authors, which 
considered gas-dynamical mechanisms that at relatively low red- 
shift are respons ible for redistrib uting previously produced met- 
als. For instance. ICora et all d2008h suggested that clumps of low- 
entropy highly enriched gas may sink in the central cluster regions, 
thereby leading to an increase of the observed emission-weighted 
metallicity. For instance, ram-pressure stripping of the interstellar 
medium (ISM) of merging galaxies has been suggested as a mech- 
anism to pollute at relative ly low redshift a metal -poor ICM with 
highly enriched gas (e.g., Domainko et al. 2006, and references 
therein), while c ausing a morpholo g ical transformation of clus 
ter galaxies (e.g., ICalura et 



lorphological transtormation ot clus- 
ai] |2007l : iRoediger & Bragged 120071) . 



Although possible evidences of r am-pressure strip ping of cluster 
galaxies have been detected (e.g. JChung et alj|2007l) the question 
remains as to whether this mechanism dominates the evolution of 
the ICM enrichment. Indeed, since ram pressure is expected to be 
more efficient in high-temperature clusters, one exp ects an increas - 
ing trend of metallicity with ICM temperature (e.g.. lRenzinill997h . 
If any, observations su ggest that hotter sy stems have a relatively 
lower metallicity (e.g.. iBaumgartner et afl l2005l) . thus suggesting 
that ram-pressure stripping is not the dominant process in enrich- 
ing the ICM. 

It is clear that understanding the history of the ICM enrich- 
ment in cosmological context, during the cluster hierarchical build 
up, requires describing in detail the gasdynamics related to the 
merging processes, while including a self-consistent treatment of 
star formation and chemical evolution. In this context, cosmologi- 
cal hydrodynamical simulations offer a unique m eans to capture in 
full detail the complexity of these processes (e.g., Valdarnini 2003 ; 
Tornatore et al.l2004lRomeo et al.l2006l : lTornatore et al.l2007l see 
Borgani et alj|2008l for a recent review). In their most advanced 
versions, chemo-dynamical simulation codes treat the production 
of different metal species, released by different stellar populations 
by resorting to detailed stellar yields, also accounting for the mass- 
dependent stellar lifetimes. 

In this paper we will present results on the ICM metal abun- 
dance from cosmological simulations of galaxy clu sters, using the 
chem o -dynamical ver sion of the GADGET-2 code dSpringel et al.l 
200 it ISpringell [2005). which has been recently presented by 
Tornat ore et al. J2007h T07 hereafter) We will compare the sim- 
ulations with observational results on the Iron abundance profiles, 
Zf , of nearby clusters, on the evolution of the ICM metallicity 
and on the Snla rates. This comparison will be performed with the 
aim of shading light on the relative role played by star formation, 
feedback processes and gas dynamics in determining the cosmic 
history of metal enrichment. 

The plan of the paper is as follows. In Section 2 we review our 
implementation of chemical evolution in the GADGET-2 code and 
present the main characteristics of the cluster simulations. Section 
3 will be devoted to the comparison between simulation results and 
observations. After comparing the profiles of the Iron abundance, 
we will concentrate on the evolution of the ICM metallicity. We 
will then compare observations and simulation predictions on the 
rate of Snla. We will draw our conclusions in Section 4. 

All values of Iron abundance that we will q uote in the follow- 
ing ar e scaled to the solar abundance value by iGrevesse & Sauvall 
dl998r) . 



2 THE SIMULATIONS 

In this letter we present a set of simulations of four massive iso- 
lated clusters, which have been identified i n a Dark-Matter onl y 
simulation having a box size 479/t _1 Mpc (Yoshid a et al1l200lh . 
performed for a flat ACDM cosmological model with Q. m = 0.3, 
hioo = 0.7, as = 0.9 and Qt = 0.04. The four extracted 
Lagrangian regions, centred on these clusters with virial massed 
in the range Af v i r =1.0-2.3 xlO 15 ft _1 M0, have been resim- 
ulated using the Z oomed Initial Condition (ZIC) technique by 
iTormen et al.l d 19971) . which allows one to increase force and mass 
resolution in the regions of interest. The high-resolution DM par- 
ticles have mass tudm = 1-13 x 10 9 /i _1 Mq, and the barionic 
particles have been added with a mass m gas — 1.7 x 10 8 ft _1 M© 
in order to reproduce the assumed cosmic barionic fraction. The 
basic characteristics of the simulated clusters are summarized in 
Tabled 

The simulations are p erformed using the hydrodynamical 
Tree-SPH code GADGET-2 dSpringelll2005h with the implemen- 
tation of chemical enrichment by T07. The Plummer-equivalent 
softening length for gravitational force is set to e — 5 h~ 1 kpc in 
physical units from z = 2 to z — 0, while at higher redshifts is 
e = 15 ft _1 kpc in comoving units. The simulatio ns include heating 
from a uniform time-dependent UV background Maardt & Madaul 
1 19961) a nd metallicity-dependent rad iative cooling based on the ta- 
bles by ISutherland & Dopitd d 19931) for an optically thin plasma. 
The process of star formation (SF hereafter) is described by th e 
sub-resolution multiphase model by Springel & Her nquistl d2003l) . 
for which the density threshold for the onset of SF is set to uh = 
0.1cm -3 . 

While the relevant features of the chemical evolution model 
are described here below, we address the reader to T07 for a more 
detailed description. Metals are produced by SnII, Snla and in- 
termediate and low-mass stars (ILMS hereafter), with only Snla 
and SnII providing energy feedback. We assume SnII to arise 
from stars having mass above 8Mq. As for the Snla, we assume 
their progenitors to be binary systems, whose total mass lies in 
the range (3-16)Mq. Metals and energy are released by stars of 
different mass by properly accounting for mass-dependent life- 
times. In this work we assum e the lifetime function proposed by 
IPadovani & Matte ucci ( 1993). We adopt the metallicity-dependent 
ste llar yields by [Wooslev & Weaver! d 1995b for SnII, the yields 
by van den Hoek & G roenewegenl d 19971) for the ILMS and by 
Thieleman n et al.l d2003l) for Snla. The version of the code used for 
the simulations presented here allowed us to follow H, He, C, N, 
O, Mg, Si and Fe. Once produced by a star particle, metals are then 
spread to the surrounding gas particles by using the B-spline kernel 
with weights computed over 64 neighbours and taken to be pro- 
portional to the volume of each particle. T07 verified with detailed 
tests that the final results on the pattern of chemical enrichment are 
rather insensitive to the weighting scheme (kernel shape and num- 
ber of neighbours) used to spread metals. 

Our simulations include the kine tic feedback model imple- 
mented by Sprin gel & Hernq uist (2003). According to this scheme, 



We define the virial mass M v i r as the mass contained within the virial 
radius R vlr . This is denned as the radius within which the average density 
p v i r is that predicted by the spherical collapse model (for the cosmology 
assumed in our simulations, p v ; r ~ 100p c , with p c the critical cosmic 
matter density). More in general, we define Ra to be the radius at which the 
density is A times the critical density p c . Physical quantities with subscript 
A are computed within R/± . 
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SnII explosions trigger galactic winds, whose mass upload rate 
is assumed to be proportional to the star formation rate, Mw = 
77 M* . Therefore, fixing the parameter r\ and the wind velocity vw 
amounts to fix the total energy carried by the winds. Our choice of 
r\ — 3 and vw = 500 k ms" 1 cor r espon ds to assume, for the initial 
mass function (IMF) bv lSalpetej dl955h . with SnII releasing 10 51 
ergs each, nearly unity efficiency in powering galactic outflows. 

In our comparison with observational data, we will first ex- 
plore the effect of changing the IMF . We u se the IMF by [Salpeter 
and that bv lArimoto & Yoshiil dl987h . for which the number 
N of stars per unit mass interval is defined as ffi(m) oc dN/Am oc 
m —( 1 + x > t w ith x = 1.35 and x = 0.95 respec tively. Further - 
more, we also use the multi-slope IMF proposed by |Kroup3 fcOOlh 
with x = —0.7,0.3 and 1.3 respectively for m ^ O.O8M0, 
0.08 ^ m < 0.5-A/q and m ^ O.5M0. Simulations based on 
the Salpeter IMF have been run for the four clusters, while only 
simulations of the g5 1 halo have been carried out for the other two 
choices of the IMF. In the following, we label the runs that use the 
Salpeter, Arimoto-Yoshii and Kroupa IMFs with 5a/, AY and Kr 
respectively. 

An important parameter entering in the model of chem- 
ical evolution is the fraction A of stars, in the mass range 
0.8-8 Mq, belonging to binary syst ems which explodes as 
Snla in the single-degenera te scenario dGre ggio & Renzini| [l983l ; 
iMatteucci & Greg gio 1986 ). For our reference runs we will use 
A = 0.1, as suggested bv lMatteucci & Gibson to re produce 

the observed ICM metallicity (see also lPortinari et alj2004l) . As we 
shall discuss in the following, the simulation with the AY IMF tends 
to overproduce Iron. In the attempt to overcome this problem, we 
also carried out a run with the AY IMF using also A = 0.05. 

Simulations of galaxy clusters, which include the scheme of 
feedback adopted here, are already known to produce an excess 
of lo w-redshift star formation, mostly a ssociated with the BCG 
(e.g jRomeo et al.|[20o3 : ISaro et alj|2006h . This recent star forma- 
tion is expected to significantly affect the history of the ICM en- 
richment. From one hand, it should provide an excess of recent 
metal production, thus possibly enhancing the enrichment at small 
cluster-centric radii. On the other hand, a recent star formation is 
also expected to lock back in the stellar phase a significant amount 
of highly enriched gas, which has shorter cooling time, thus leav- 
ing in the hot ICM only relatively metal-poorer gas. In order to 
quantify the effect of recent star formation of the ICM enrichment 
history, we have also simulated the 5a/ version of the g5 1 cluster 
by switching off radiative cooling and star formation below 2 = 1, 
considering both the case in which already formed stars keep pro- 
ducing metals with the appropriate lifetimes (CS run) and the case 
in which also the metal production is stopped at the same redshift 
(CMS run). While this prescription of suppressing low-redshift star 
formation and metal production is admittedly oversimplified, it al- 
lows us to address the following questions: (i) to what extent the SF 
excess in simulations affects the enrichment evolution of the ICM? 
( ii) which is the role of gas-dynamical processes in redistributing at 
relatively low redshift the metals that have been produced at earlier 
epochs? 



3 RESULTS AND DISCUSSION 

3.1 Metallicity profiles of nearby clusters 

The radial profiles of the metal abundance provide a very impor- 
tant record of the chemical enrichment process in galaxy clusters. 
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Table 1. Characteristics of the simulated clusters at z = 0. Column 1: 
cluster name; column 2: virial mass (units of 10 15 /i _1 M©); column 3: 
virial radius (units of fe~ 1 Mpc); column 4: spectroscopic-like temperature 
within i? v i r (keV, see Mazzotta et al. (2004), for its definition). 



Indeed, they are determined by the distribution of cluster galax- 
ies, where most of the metals are produced, by the mechanisms re- 
sponsible for their transport and diffusion from the star-forming re- 
gions (i.e., galactic ejecta, ram-pressure and viscous stripping, etc.) 
and by other gas-dynamical processes which redistribute them on 
larger scales (e.g. turbulence and sinking of enriched low-entropy 
gas). Here we compare the profiles of the Iron abundance of simu- 
lated galaxy clusters at 2 = with the observational results from 
Chandra data of a sam ple of nearby relaxed clusters analysed by 
IVikhlininetaT] ( l2005h . 

In FigureQ]we compare the profiles of Zf c from our simulated 
clusters with the observed profiles of 8 clusters h aving temperature 
above 3 keV (see Table 1 in Vikhl inin et ai]2 005). We point out that 
the analysis o f lVikhlinin et alJ ((2005) provided information on the 
total ICM metallicity, i.e. without distinguishing the contribution 
from different chemical species. However, at the typical tempera- 
tures of these clusters and for the typical energy r ange where the 
spect ral analysis was performed (0.6-10 keV; see I Vikhlinin et all 
2005), this observed metallicity is largely dominated by Iron. We 
want to stress the fact that the simulated clusters are dynamically 
relaxed (with the last major merger undergone before z = 0.5) 
and therefore suitable for the comparison with this set of observed 
clusters. 

As shown in the left panel of FigQ] the simulations based on 
a ISalpeterl dl955l) IMF produce profiles which are in reasonable 
agreement with observations. The scatter among the four simulated 
clusters is quite small, with some increase in the central regions, 
7? ^ 0.1 i^soo- Although observations seem to have a larger scatter, 
it is not clear how much observational uncertainties contribute to it. 
Changing the IMF (right panel of Fig. [T} clearly turns into a change 
of the overall amount of the Iron abundance at all radii, with both 
the Kroupa and the Ay IMFs producing too high profiles. The larger 
amount of Iron found for these two IMFs is due to the fact that, 
once normalized, they both predict a larger number of supernovae 
contributing to the Iron production, with respect to the Salpeter one. 

Besides producing more SnII, the AY and Kr IMFs also pro- 
duce a larger number of Snla, since there is a significant overlap 
between the mass range relevant for Snla and the mass range where 
these two IMFs are higher than the Salpeter one. 

As for the relative roles of Snla, SnII and ILMS in the ICM 
enrichment, we verified that Snla contribute for about 70 per cent 
of the Iron contained in the diffuse medium within R500 for the 
Salpeter IMF. This fraction decreases to about 65 per cent for the 
Kroupa IMF and to about 55 per cent for the Arimoto-Yoshii IMF. 
Since Snla provide a major contribution to the Iron production, our 
results are quite sensitive to the choice of the fraction A of stars 
in binary systems. As a matter of fact, this fraction can be consid- 
ered as a free parameter in a model of chemical evolution. Follow- 
ing a phenomenological approach, for each choice of the IMF its 
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Figur e 1. The Iron a bundance profiles of simulated clusters at z = compared to the Chandra observed profiles for 8 nearby clust ers with T > 3 keV 
I Vikhli niiTet alj2005l) (open circles with errorbars). In the left panel the solid curves correspo nd to the four clusters simulat ed using thelSatoetej jl955 ) IMF. 
The right panel compares the results of the g51 cluster simulated using the IMF by Salpetei ( 1955) ( solid) , bv lArimoto & Yo shii (1987) with two different 
values for the fraction of binary stars (A = 0.1: dot-dashed; A = 0.05: short dashed) and by Kroupa 12001) (long dashed). 



value is determined by the requirement of reproducing some obser- 
vational data. In our case, we note that decreasing A from 0.1 to 
0.05 induces a significant decrease of the Zp c profile. This sort of 
degeneracy between the IMF shape and the fraction of binary stars 
can be broken by looking at the relative abundance of a elements 
with respect to Iron. For instance, since Oxygen is essentially pro- 
duced by SnII, we expect a top-heavier IMF to provide values of 
O/Fe higher than for a top-lighter IMF. If we suppress the num- 
ber of Snla for the top-heavier IMF, by decreasing the value of A, 
we further increase the O/Fe ratio, thus allowing to distinguish this 
case from that of a top-lighter IMF with a higher A. We deserve 
a forthcoming paper to a detailed comparison of simulations with 
observations of relative abundances for nearby clusters. 

The runs with the Salpeter IMF provide results in closer agree- 
ment with the Chandra data, although in all cases the profiles of the 
simulated clusters are somewhat steeper than the observed ones, 
with negative gradients extending at least out -R5 00 and beyond. 
This r esult is at variance with the recent claim by ISnowden et al.l 
d2008h who found no evidence for the presence of metallicity gra- 
dients at scales ^ O.I-R500 from the analysis of a catalog of 70 
clusters observed with XMM-Newton. In the next section we will 
compare simulated and observed results on the evolution of the 
ICM metallicity at small radii, where the simulated and the ob- 
served metallicity gradients are in reasonable agreement. If con- 
firmed by independent analyses, the lack of abundance gradients 
at relatively large radii will provide a non-trivial constraint for 
chemo-dynamical models of the ICM enrichment. 

Limited numerical resolution could lead to an underestimate 
of high-redshift enrichment from a pristine population of relatively 
small under-resolved galaxies. This high-2 enrichment should be 
rather uniform and, therefore, should soften the metallicity gradi- 
ents. Indeed, T07 found that increasing resolution provides pro- 



gressively shallower metallicity profiles. However, the effect is vis- 
ible only at radii C 0.5-Rsoo , while being negligible at smaller radii, 
which are dominated by the star formation associated to the BCG. 
Another possibility to soften met allicity profiles can be pr ovided 
by AGN feedback. For instance, [Bhattacharva et all {2007) anal- 
ysed cosmological simulations of galaxy groups, which include the 
effect of energy feedback from gas accretion onto black holes. They 
found that the effect of this feedback is to redistribute the hot gas, 
driving it from the inner regions, where it should be more enriched, 
to the outer part of the halo, and to lower the star formation in 
the inner region. Siiac ki & Springell J2006ah used a similar feed- 
back scheme, in which energy is used to trigger the formation of 
high entropy bubbles. These bubbles rise buoyantly in the ICM, 
giving ris e to a redistribution of the central metal-enriched gas 
(see also Roed igeretZi r2007). Clearly, in this case the request 
is that the redistribution of metals should not be so efficient as to 
destroy the metallicity gradients in the central cluster regions (e.g., 
iBohringer eTal]|2004l) . 



3.2 Evolution of the ICM metallicity 

In this Section we compare the simulation predictions on the evo- 
lution of the ICM m etallicity with the observational results by 
iBalestra eTa l. (2007). These authors analysed Chandra observa- 
tions of 56 clusters at z > 0.3 (with the addition of XMM-Newton 
observations for clusters at z > 1) having temperatures above 3 
keV. They measured the metallicity in the central regions, with a 
typical extraction radius of 0.15-0.3 Jiiso, chosen object-by-object 
so as to maximize the signal-to-noise ratio (see also lMaughan et all 
20081 r or a similar an alysis). For the low-redshift reference value, 
Balestra et al. I J2007I) combined this set of distant clusters with a 
mix of cool-core and non cool-core clusters at lower redshift. They 
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Figure 2. The comparison between observations and simulations for the evolution of the Iron abundance, Zp a . Observational results from [Balestra et alj 
1 2007) are shown with open circles, with errorbars corresponding to the la uncertainty in the combined spectral fit performed for all the clusters falling within 
each redshift bin. The shaded area is the r.m.s. scatter among the measured metallicities within the same redshift intervals. Left pa nel: the dependence of the 
simulation results on the stellar IMF. The filled squares show the average over the simulated clusters, assuming a Salpeter 11955) IMF (Sal), with errorbars 
indicating the r.m.s. scatter over the four objects. For the g51 cluster only, the open squares are for the run with Sal petea g 19550 IMF, the filled and open 
triangles are for the Arimoto & Yoshiil J 1 9871) IMF (AY) with A = 0.1 and A = 0.05 for the binary fractions, respectively, while the filled circles are for the 
run with the Kroupa (2001) IMF (Kr). Right panel: the effect of stopping star formation and metal production at low redshift on the evolution of Zp c for the 
g51 cluster, using the lSalpeteil Jl955T > IMF. Results for the reference (Sal) run are shown with the open squares. The filled circles are for the simulation with 
radiative cooling and star formation stopped at z = 1 (CS), while the filled triangles are for the run in which also the metal production is turned off at z = 1 
(CMS). 



also pointed out that the observed decrease of Zf c with redshift is 
not induced by a decrease of the fraction of cool-core clusters in 
the past. Therefore, we expect that no significant bias is introduced 
when comparing the observed evolution with that traced by our set 
of relaxed simulated clusters. Since it is quite difficult to define a 
common extraction radius for the observed clusters, we decided to 
adopt a value of 0.27?iso in the analysis of the simulated clusters. 
We verified that our conclusions are left unchanged by varying this 
radius in the range 0.15-0.3 Riso- 

In the left panel of Figur e [2] we compare the observational 
results by Bales tra et al.1 d2007t) with the predictions of our simu- 
lations for different choices of the IMF. Observations and simula- 
tions are compared here by using the emission-weighted definition 
of metallicity, with emissivity of each gas particle computed in the 
0.5-10 keV energy band. In principle, this comparison would re- 
quire extracting synthetic spectra from the simulated clusters and 
then measure the metallicity by fitting these spectra to a single- 
temperature and single-metallicity plasma models, as done in the 
analysis of observatio nal data. An anal ysis of this type has been 
recently presented by Rasi a et all J2007h and showed that, at least 
for Iron, the emission-weighted estimator gives results quite close 
(within about 10 per cent) to those obtained from the spectral- 
fitting analysis. 

Interestingly, we note in all simulations a significant increase 
of metallicity in the cluster central regions below redshift unity. 
The runs based on the 5a/ IMF provide results close to observa- 
tions. On the contrary both the Kr and the AY IMFs predict too high 



abundances for the g51 cluster at all redshifts, with a very strong 
evolution at z^, 0.5 (we note that g51 is the cluster with the highest 
metallicity at z = 0.25 and 0.5, among the four simulated objects). 
This higher abundance for the AY and Kr IMFs is in line with the 
correspondingly higher profiles found at z = (see Fig.[T). Again, 
decreasing the binary fraction in the AY run to A = 0.05 causes a 
significant decrease of the Iron abundance at all redshifts. 

At lower redshifts, z > 0.5, g51 shows a slower evolution, 
even followed by an inversion below z ~ 0.25. The reason for this 
inversion lies in the quiet accretion history of this clusters below 
z — 0.5. Since no highly enriched gas clumps, associated to major 
merger event, reached the central regions of g51 since z — 0.5, the 
only gas accreted there is characterized by a relatively low metal 
abundance. For instance, in the Sal run of g51 we verified that 
beetween z = 0.25 and z — about 6 x 10 12 M© of gas, hav- 
ing an average metallicity value Zf& ~ 0.64 Zfc,q, flowed out of 
0.2i?i8o- In the same redshift interval, about 1O 13 M0 of gas was 
accreted, with an average metallicity Zfe — 0.42 Zfe,q- There- 
fore, in the absence of large metal-enriched clumps reaching the 
cluster centre, gas mixing leads to the accretion of relatively metal- 
poorer gas, thus turning into a decrease of Zp c . 

Although the agreement between the runs based on the 
Salpeter] d 1955b IMF and observations is rather encouraging, the 
question remains as to whether the positive evolution seen in the 
simulations is just the spurious product of the excess of star forma- 
tion taking place in the central cluster regions. In order to address 
this question, we compare in the right panel of Fig|2]the evolution 
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Figure 3. Projected maps of the emission-weighted Iron abundance for simulations of the g5 1 cluster, using the IMF by Salpeter] ll955l) . for the reference run 
(upper left), for the run with cooling and star formation stopped at z = 1 (CS run; upper right) and for the run in which also the metal production is stopped 
at z = 1 (CMS run; lower left). Each map is 2i? v i r a side. 



of Zyc for g51 when stopping star formation and/or metal produc- 
tion at low redshift. Quite remarkably, halting star formation below 
z = 1 while allowing already formed stars to keep releasing metals 
(CS run) has the effect of strongly increasing the positive evolution 
of Zyc in the central region of g51, which turns out to be over- 
enriched by z — 0. This leads to the counter-intuitive conclusion 
that the lack of low-z star formation should generate an increase of 
the enrichment of the hot gas. In order to investigate the origin of 
this increase, we show in Figure[3]the emission-weighted metallic- 
ity maps of the reference Sal run of g5 1 at z — 0, along with those 
of the CS and CMS runs. Quite apparently, the metal distribution in 
the CS simulation is more clumpy than in the reference run. At the 
cluster centre, a high Zf c is clearly visible, which boosts the cen- 
tral emission weighted metallicity shown in Fig[2] Indeed, while 
the emission-weighted Zf c increases by about a factor two within 
0.2i?i8o, we verified that the mass-weighted estimate within the 
same radius increases only by about 10 per cent. In the reference 
run, the metals released in the high density clumps disappear from 
the hot diffuse medium due to the efficient gas cooling. As a re- 
sult, the reference run has a globally higher level of diffuse enrich- 
ment, but a lower level of enrichment inside the high-density gas 
clumps, which dominate the emission-weighted estimate of Zf c . 
These results demonstrate that the strongly positive evolution of 



the emission-weighted metallicity in the CMS run is driven by the 
accretion of highly enriched dense clumps. 

Inhibiting also the production of metals below redshift unity 
(CMS run) allows us to characterize the role played by gas- 
dynamical processes in redistributing metals produced at higher 
redshift. As shown in the bottom-left panel of Fig[3] metal clumps 
are less pronounced than in the CS run. The global enrichment level 
of the ICM is now significantly lower than in the reference run, al- 
though an enhancement in the innermost regions is still visible. The 
resulting mass-weighted metallicity within 0.2i?i8o at z = de- 
creases by ~ 60 per cent with respect to the reference run. There- 
fore, the stability of the emission-weighted metallicity is due to the 
competing effects of a more clumpy distribution of metals and of a 
decrease of the overall ICM metal budget. 

The maps of FigfJ] also illustrate the role of gas-dynamical 
effects in redistributing highly enriched gas. Merging clumps 
within the cluster virial region leave behind them over-enriched 
tails of stripped gas, which is tempting to explain as due to 
ram-pressure stripping. However, a significant contribution could 
well be provided by viscous stripping. Since the SPH scheme is 
known to be generally characterized by a large numerical viscos- 
ity, this may induce an exc ess of gas stripping from merging halos. 
ISiiacki & Springell {2006b ) showed that the effect of including the 
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Spitzer-Braginskii viscosity in the SPH, on the top of the numeri- 
cal viscosity, is indeed that of furthe r increasing gas stri pping from 
merging halos. On the other hand, iDolag et alJ J2005t) discussed 
an SPH scheme of reduced viscosity. In this case, the increase of 
the "turbulent" stochastic gas motions should provi de a more effi- 
cient diffusion of metals from star-forming regions (Rebus co et al.l 
120051) . while making viscous stripping less efficient. Although it is 
beyond the aim of this paper to carry out an accurate analysis of the 
effect of viscosity on the pattern of the ICM enrichment, there is no 
doubt that this aspect deserve an accurate in-depth investigation. 



3.3 The Snla rate 

The supernova rate represents a useful diagnostic to link the ob- 
served evolution of the ICM metallicity to the past history of star 
formation and to shed light on the relative contribution of Snla 
and SnII in releasing metals. In particular the ratio between the 
Sn rate and the B-band luminosity, the so-called SnUs, can be 
used to distinguish the relative contribution of Snla, which form in 
binary systems of stars with masses in the range (0.8-8) Mq and 
the short-living massive stars that contribute substantially to the 
B-band luminosity of galaxies. In this section we present a com- 
parison between the results of our simulate d clusters and observa - 
tional data of SnUfl in gal axy clusters fr om Gal-Yam et alJ d2002l) . 
iMannucci et"ai1 d2008h and lSharon etal] d2007l) . 

The simulation analysis finalized to compute the SnUs pro- 
ceeds as follows. For each star particle we know its formation red- 
shift and metallicity. Given the IMF and the lifetime function, this 
allows us to compute the rate of Snla exploding in each such par- 
ticle. Furthermore, using the spectrophotometric GAL AXE V code 
( iBruzual & Charlotll2003h we also compute the B-band luminosity 
of each star particle, which is treated as a Single Stellar Popula- 
tion (SSP). Once Snla rates and luminosities are computed for all 
the star partic les, we run the SKID substructure-finding algorithm 
dStadelll200ll) on their distribution to identify galaxies as gravi- 
tationally bound groups of stars. All the star particles not bound 
to ga laxies take part of th e intra-cluste r diffuse stellar c omponent 
(e.g jMurante et alj|2007h . We refer to lSaro et al.l d2006l) for a de- 
tailed description of the procedure to identify galaxies and assign 
broad-band luminosities to them. In order to reproduce the obser- 
vational procedure, we compute the SnUs values by also including 
the contribution of the Snla arising from diffuse stars, while the B- 
band luminosity is computed by including only the contribution of 
the identified galaxies. 

In the left panel of Figure|4]we compare the SnUs values from 
the simulations with different IMFs with observational data. In per- 
forming this comparison one potential ambiguity arises from the 
definition of the extraction radius, within which luminosities and 
Snla rates are measured in observations, since different authors use 
different aperture radii. To address this issue we computed SnUs 
in the simulations within R v i r and verified that the results are left 
unchanged when using instead R500 . 

Observational data show a declining trend at low redshift. This 
is generally interpreted as due to the quenching of recent star for- 
mation, which causes the number of Snla per unit B-band lumi- 
nosity to decrease after the typical lifetime of the Snla progenitor 
has elapsed. On the other hand, our simulations predict a rather flat 
evolution of the SnUs, independently of the choice for the IMF, 
which is the consequence of the excess of low-redshift star for- 
mation. The runs based on the Salpeter and on the Kroupa IMF 
produce very similar results. Although the Kroupa IMF produces 
a higher rate of Snla, due to its higher amplitude in the (1-8) Mq 



stellar mass range, this is compensated by the higher values of Lb ■ 
These two IMFs both agree with the observational data at z si 0.3 
within the large observational uncertainties, while they overpredict 
the rates measured for local clusters. Although the excess of recent 
star formation in the c entral regions of our simulated clusters pro- 
duces too blue BCGs 1 Sam et al. 2006;), the large number of Snla 
associated to this star-formation overcompensate the excess of blue 
light. 

As for the simulation with the Arimoto-Yoshii IMF, it predicts 
an even higher SnUs at low redshift. As shown in the right panel 
of FigO decreasing the binary fraction to A = 0.05 decreases 
the value of the SnUs by more than a factor 2. While this helps 
in reconciling the simulation results with the low-redshift data, it 
introduces a tension with the data at z ~ 1. 

Truncating star formation at z — 1 (right panel of FigO has 
the desired effect of decreasing the value of SnUs below z = 0.5. 
Quite interesting, for 0.5Js z^, 1 the decreasing trend of the Snla 
rate is compensated by the corresponding decrease of the B-band 
luminosity, while it is only at z > 0.5 that the decrease of the Snla 
rate takes over causes the decrease of the SnUs values. 



4 CONCLUSIONS 

We have presented results from cosmological SPH hydrodynam- 
ical simulations of galaxy clusters with the purpose of charac- 
terising the evolution of the chemical enrichment of the intra- 
cluster medium (ICM) out to redshift z ~ 1. The simulatio ns have 
been performed with a version of the GADGET-2 code (jSpringel 
12005b. which includes a detailed model of chemical evolution 
dTornatore et alj|2007l) . Our simulations have been performed with 
the purpose of investigating the effect of changing the chemical 
evolution model and the effect of suppressing star-formation at 
z < 1. The main results of our analysis can be summarized as 
follows. 

(a) T he Iron abundan ce profiles provided by simulations based 
on the ISalpeted dl955l) IMF are in reasonable agreement with 
the results from Ch andra observations of nearby clusters by 
IVikhlinin et alJ d2005l) a t 7?< 0.2 R 500 . Simulations based o n the 
IMFs bv lKroupal d200lh (Kr) and lArimoto & Yoshiil dl987l) (AY) 
both predict too high a normalization for these profiles. However, 
reducing the fraction of stars assumed to belong to binary systems 
suppresses the enrichment level, thus alleviating the disagreement 
of a top-heavy IMF with the observed Zf c profiles. Our simula- 
tions always predict negative metallicity gradients extending out to 
-R500 and beyond, possibly in disagreement with XMM-Newton 
measurements of the Iron metal abundance at relatively large radii 
dSnowden et"ai]|2008h . 

(b) All our simulations predict a positive evoluti on of the cen- 
tral Iro n abunda nce, comparable to th at observed by Bales tra et al.l 
d2007l) (see also lMaughan et alj |2008). Using a Salpeter IMF also 
provide an enrichment consistent with observations, while the 
Kroupa and Arimoto-Yoshii (AY) IMFs overpredict the enrichment 
level at all redshifts. Again, this disagreement can be alleviated by 
decreasing the fraction of binary systems. It is worth reminding that 
the observed evolution of the Iron abundance is traced by using a 
mix of cool-core and non cool-core clusters, while our simulated 
clusters are all dynamically relaxed. Clearly, a fully self-consistent 
comparison would require simulating a representative population 
of clusters, having a variety of morphologies and dynamical states. 

(c) Stopping cooling and star formation at z — 1 (CS run) has the 
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panel: the effect of changing the IMF. The open squares ar e for the Sal peteJ ll955t) IMF, the open triangles are for the top-heavy IMF by Arim oto & Yoshiil 
1 19871) and the open circles for the IMF by iKroupa 1200 lh - For the Salpeter IMF, the shaded area show the r.m.s. scatter evaluated over the four simulated 
clusters, while for the other two IMFs only the result for the g5 1 cluster is shown. Right panel: the effect of suppressing low-redshift star formation and of 
changing the binary fraction on the SnUs evolution of the g51 cluster. The open squares and the open circles are for the reference run with lSalpeteJ 1 1955) 
IMF and for the same ru n with cooling and star formation stopped at z = 1 (CS run), respectively. The filled and the open triangles are for the runs with 
lArimoto & Yoshiil 1 19871) IMF, using A = 0.1 and A = 0.05 for the fraction of binary stars, respectively. 



effect of producing a too strong positive evolution of the emission- 
weighted metallicity. Indeed, in the absence of star formation all the 
metals released at z < 1 by long-living stars are no longer locked 
back in the stellar phase. As a result, metallicity is enhanced inside 
high-density halos and in the central cluster region. The clumpy 
metal distribution boosts the emission-weighted abundance esti- 
mate. This leads to the somewhat counter-intuitive conclusion that 
suppressing recent star formation has the effect of enhancing the 
positive evolution of the ICM metallicity. 

(d) A comparison of the Snla rate per unit B-band luminosity, 
SnUs, show that our simulations are generally not able to repro- 
duce the observed declining trend at low redshift. This result is ex- 
plained by the excess of recent star formation taking place in the 
central regions of galaxy clusters. Indeed, excising star formation 
at z < 1 produces an evolution of SnUs which is consistent with 
the observed one. 

Cluster simulations which only include stellar feedback, like 
those presented here, are well known to be at variance with a num- 
ber of observations, such as the temperature profiles in the cool 
core regions and an large excess of recent star formation in the 
BCG. Our prescription to quench recent star formation is admit- 
tedly oversimplified. A more realistic treatment would require in- 
troducing energy feedback from gas accretion onto super-massive 
black holes, which sel f-consistently follow the h ierarchical build- 
up of the cluster (e.g.. lSiiacki & Springel 2006a). Still, our results 
highlight that the positive evolution of the metal abundance in the 
central regions of simulated clusters can not be simply interpreted 
as a consequence of an excess of low-redshift star formation. In 



fact, the evolution of the metallicity pattern is driven by the com- 
bined action of the gas-dynamical processes, which redistribute 
already enriched gas, and of star formation, which acts both as a 
source and as a sink of metals. While hydrodynamical simulations 
probably provide the most complete interpretative framework for 
observations of the history of the ICM enrichment, they have still 
to improve in the numerical accuracy for the description of relevant 
physical processes. There is no doubt that the ever increasing super- 
computing power and efficiency of simulation codes should be par- 
alleled by a comparable advance in the reliability of the numerical 
description of the relevant astrophysical and gas-dynamical pro- 
cesses. Only in this way, simulations will become the standard tool 
to link ICM observations to the global picture of cosmic structure 
formation. 
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